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I. INTRODUCTION 

The spin-orbit potential for a nucleon moving in a nu- 
cleus plays a decisive role in the nuclear shell structure. 
Though the phenomenological aspect of the spin-orbit 
potential has been clarified rather well for stable nuclei, a 
full account of its origin is still under active study. Impor- 
tant information on the spin-orbit potential is obtained 
through the analysis of elastic scattering observables, es- 
pecially a vector analyzing power in the elastic scattering 
of polarized protons from the nucleus. 

Recent developments of experimental techniques have 
provided us with a polarized proton target and have made 
it possible to measure not only the differential cross sec- 
tion but also the vector analyzing power for the elastic 
scattering of polarized protons from unstable nuclei in the 
inverse kinematics. The data have recently been taken 
for 6 Hc [l| and 8 He at 71 MeV/nucleon. A good ex- 
ample of a two- neutron halo nucleus is 6 He, where the 
neutron density extends far out in distance. This unique 
feature of its structure is expected to show up in the vec- 
tor analyzing power because the spin-orbit potential is 
primarily sensitive to the surface of the nucleus. Though 
less pronounced than 6 He, 8 He has an extended neutron 
cloud as well and its clastic scattering observables are 
interesting. 

There have been so far a few theoretical studies on the 
optical potential and observables for p+ 6 ' 8 He elastic scat- 
tering. The single scattering approximation to the mul- 
tiple scattering expansion was employed in the p+ 8 He 
case Q . The predicted angular distribution of the vector 
analyzing power shows a peak at about # c .m.=46°, dif- 
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ferent from the data Q. The p+ 8 He elastic scattering 
angular distribution was analyzed in the eikonal model 
to examine the sensitivity to the matter distribution of 
8 He [11 ■ The p+ 6,8 He elastic scattering observables were 
calculated in a full- folding optical model Q. The calcu- 
lated vector analyzing powers do not agree with exper- 
iment P, [|| • Very recently the neutron pickup coupling 
with the elastic channel has been studied to see its effect 
on proton scattering from 6 He @ . As noted above, none 
of the calculations succeeds in reproducing the vector an- 
alyzing powers for p+ 6,8 Hc scattering. It is reported that 
a search for the optical potential parameters of the spin- 
orbit part leads to a shallow and long-ranged spin-orbit 
potential for 6 He [l|. 

The purpose of this paper is to analyze the elastic scat- 
tering observables, the differential cross section, the vec- 
tor analyzing power and the spin-rotation function, for 
protons scattered from the He isotopes including 4 He Q ■ 
The central part of the optical potential for the proton 
is calculated in the framework of Glauber or the eikonal 
model 0. The inputs needed in the calculation include 
only the ground state wave functions of the He isotopes 
and the nucleon-nucleon scattering amplitude, or more 
precisely its Fourier transform, the nucleon-nucleon pro- 
file function. The spin-orbit potential is constructed us- 
ing a derivative of the central part of the optical poten- 
tial. There are at least three noticeable advantages of 
the present approach: Firstly it is logically very simple, 
and nevertheless contains nucleon-nucleon multiple scat- 
terings to all orders. Secondly the wave function of the 
projectile nucleus itself can be employed, though in an 
approximate version of the present approach the projec- 
tile density may be used instead of the wave function. 
Thirdly the optical potential obtained takes account of 
breakup effects of the projectile without recourse to labo- 
rious calculations with continuum discretization P, [loj , 
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which enables us to discuss easily a dynamic polarization 
potential (DPP). The breakup effect is a vital ingredient 
that should be taken care of for the optical potential of 
a weakly bound nucleus like 6 He. 

In Sec. |TI]we present a formulation needed in our ap- 
proach. We show in Sec. Ill Al how to construct the opti- 
cal potential, its approximate version and a relationship 
between those potentials and a folding potential. We de- 
fine our spin-orbit potential in Sec. IIIBI together with 
the elastic scattering observables. In Sec. Ill Cl we outline 
the variational Monte Carlo method that is employed to 
generate the ground state wave functions of the He iso- 
topes. Results of calculation are presented in Sec. IIIII 
First we compare the elastic scattering observables at 
71 MeV/nucleon with experiment in Sec. IIII Al The an- 
gular distributions of the differential cross section and 
vector analyzing power for p+ 4 Hc scattering at higher 
energies are compared to available data in Sec. IIII Bl Pre- 
dictions for p+ 6,8 He scattering at intermediate energy are 
also given. Conclusions are drawn in Sec. IIVI 

II. FORMULATION 
A. Optical potential in the eikonal approximation 

Let us assume that the projectile nucleus moves in the 
z direction and impinges on a proton target with a ve- 
locity v. Let R stand for the relative distance vector be- 
tween the projectile and the proton. Under the eikonal 
approximation the x, y component of R, denoted b, turns 
out to be just a parameter, the impact parameter. The 
interaction v p n between the ith nucleon of the projectile 
and the proton gives rise to a multiplicative phase fac- 
tor e lXpJV ( b-Si ' [8| that modifies the wave function of the 
projectile, where Si is the x,y component of the position 
vector Ti of the ith nucleon relative to the center of mass 
of the projectile. The phase XpN is related to v p n by 



X P N{b) = -j-j v pN (Vb 2 + z 2 )dz. (1) 

Each nucleon contributes the position-dependent multi- 
plicative phase factor. A key quantity to describe the 
proton-projectile elastic scattering is given by the eikonal 
phase 

e *»(&) = ( $o | e iE(6, Sl> ..., 8A )[ $o ^ 

where E(b, Si, . . . , sa)= Y,t=iX P N(b - s ?; ) is the total 
phase, and <£>o is the ground state wave function of the 
projectile nucleus. 

As is well-known 0, [lj| , the optical phase shift func- 
tion ^ obtained in the eikonal approximation includes 
the effects of coupling with excited states or breakup con- 
tinuum states. To make this point clear, we define the 
average of the total phase 

Xp(b) = ($ |B(Mi,...,« A )|#o>. (3) 



Hereafter the projectile nucleus is assumed to be spheri- 
cal, so that both Xe(&) and xf(£>) become a function of 
b = \b\. Using Eqs. ((T|) and we find that xf(&) is the 
phase shift function corresponding to the single folding 
potential Uf(R)= J pN(r)v P N(\R — r\)dr 

i r + °° i 

Xf(&) = / U { (Vb 2 + z 2 )dz, (4) 

where pw(r) is the nucleon density of the projectile nu- 
cleus. The eikonal phase Xe(^) is expressed as follows 

XE(b) = Xf(&) - iln(* | c^ s --^)-»( fc »|$ ). (5) 
According to Glauber ||, a potential U C (R) 

U C (R) = —--£- [°° XE (VBF+^)dx (6) 
7r R dR J 

produces the phase shift function xe(^) in the eikonal 
approximation. The potential U C (R) differs from Uf(R) 
in that the underlying phase shift function of the for- 
mer contains the second term of the right-hand side of 
Eq. (j5|. The term can be discussed in a cumulant expan- 
sion 0, E3] that involves the fluctuation of the higher or- 
der cumulants ($ | (2(6, Si, . . . , sa) - xf(&))"|^o)- The 
difference between U C (R) and the folding potential Ui(R) 
is the DPP. As is clear from the above derivation, the 
DPP is evaluated in the eikonal approximation without 
an explicit invoking of the couplings with excited and 
continuum states. 

It should be stressed that the central part of the op- 
tical potential U C (R) can be obtained in a unified way 
independently of the projectile nucleus. 

The above formulation has successfully been applied 
to a study of the breakup effects of weakly bound pro- 
jectile nuclei, e.g., d scattered from 58 Ni Q and 6 He 
scattered from 12 C [ll|. In these applications the tar- 
get is not a proton but a composite nucleus. It is 
treated as an absorbing point particle and then it is 
possible to apply exactly the same formulation as above 
by adopting appropriate nucleon-targct optical potentials 
for v p n- The quality of such calculations is tested by 
comparing to other calculations that explicitly include 
the breakup channels in continuum discretized coupled- 
channels method [H], [l3| . 

For the proton target, v p n stands for the proton- 
nucleon potential. Any operator dependence of the po- 
tential has to be avoided because otherwise the eval- 
uation of Eq. ([2]) together with Eq. (TTJ) is impossible. 
Instead of looking for some effective forces that have 
no operator dependence, we follow a simple procedure 
here. We introduce the proton-nucleon profile function 
T pN (b), which is equal to 1 -c 4 ^"^. With T pN {b), the 
proton-nucleon scattering amplitude f P N{9) is given in 
the eikonal approximation by 

f P N(e) = ^ Je-^ b r pN (b)db, (7) 
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where K is the wave number and q is the momen- 
tum transfer, g=2i4Tsin|, and r p jy(b) is conveniently 
parametrized as 



r P iv(b) 



1 - ia pN 
4tt(3 pN 



> p n 1 



b 2 

2 PpN 



(8) 



where a P Ar is the ratio of the real to the imaginary part 
of the pN scattering amplitude in the forward direction, 
f5 p N is the slope parameter of the pN elastic differential 
cross section, and a p °^ is the pN total cross section due to 
the nuclear pN interaction. We use the parameter values 
tabulated in Ref. [l4|. The difference between pp and 
pn interactions is taken into account in what follows by 
extending Eq. © to 



r P Ar(fe) = 8jy p r p p(b) + 8jy n r pn (b). 



(9) 



In this way we bypass the direct use of the nuclear force 
in calculating xe(^)- 

The eikonal phase Xe(&) and the average total phase 
Xf(^) are expressed in terms of the profile function as 
follows 

A 

X E (6) = -ilog($ | n C 1 -!>(*>- I$o}, (10) 

i=l 

Xp(6) = -i f PN(r) log (l - T pN (b - s)) dr. (11) 

Both neutron and proton densities are employed in cal- 
culating xf(&)- As will be explained in Sec. Ill CI xe(&) 
is obtained with a Monte Carlo integration. It turns out 
that the present approach together with the Monte Carlo 
integration is very versatile for calculating the eikonal 
phase and the corresponding optical potential. A simpler 
calculation is to take the leading order of the cumulant 
expansion, leading to 



PN(r)T pN (b - s)dr. 



(12) 



This provides us with a reasonable approximation to the 
full phase for the nucleon-nucleus scattering, and is of- 
ten adopted for the evaluation of the phase shift function 

I 



for the scattering between complex nuclei [14| - |16l ] . Equa- 
tion © is used to obtain those potentials which generate 

the phases xf(&) and Xe an d they are compared with 
U C (R). 



B. Elastic scattering observables 

The potential U C (R) constructed in this way is central 
and includes the effect of breakup of the projectile in 
the eikonal approximation. The spin-orbit term of the 
optical potential is introduced as 



u so (R) = v so xl^u c (R), 



(13) 

where A w is the pion Compton wavelength. The real and 
imaginary parts of the spin-orbit potential are thus ob- 
tained from the real and imaginary parts of U c , respec- 
tively and the constant V^ that determines the spin-orbit 
strength is allowed to be different depending on the real 
and imaginary parts of the spin-orbit potential. 

Proton-nucleus scattering is described in a partial wave 
expansion with the potential U C (R) + I ■ crU so (R). The 
scattering amplitude becomes an operator in the spin 
space 



f(0) = f(0) + ig(6)<r ■ h 



(14) 



with a unit vector n perpendicular to the scattering plane 



kxk' 
Ifc x fcV 



(15) 



where k and k' are the momenta in the center of mass 
before and after the scattering. Denoting the S matrix 
by 5^ = exp(2ii5j ± ) where Sf 1 are complex phase shifts for 
the potential U C (R) + lU so (R) or U C (R) - {I + l)U BO (R), 
respectively, we obtain 



f(9) = fc(0) + 2ijfe£[(*+ ^ J ) + l ( S f !)] e2lS?Pli - c 
i 



(16) 
(17) 



where fc(&) is the Coulomb scattering amplitude, 6f is 
the Coulomb phase shift, and PHO) is the associated Leg- 
endre polynomial. A deviation of the Coulomb potential 
from that of a point charge is included in the calculation 
but its effect is very small. 



For a fixed angle 9 the scattering amplitude / is deter- 
mined by three real quantities. They can be convention- 
ally chosen as the differential cross section dcr/dfl, the 
vector analyzing power A y , and the spin-rotation func- 
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tion Q: 



The two-body correlation operator has the structure 



% = \m? + \m?> 



A y {6) 
Q{0) 



2Rc(/(%*(g)) 

\f(0)\ 2 + \g(0)\ 2 

2Im(/(fl)g*(0)) 

\f(0)\ 2 + \g(O)\ 2 



(18) 
(19) 
(20) 



C. Variational Monte Carlo wave function 

The wave functions of He isotopes used in this work are 
taken from variational Monte Carlo (VMC) calculations 
for a Hamiltonian consisting of nonrelativistic nucleon ki- 
netic energy, the Argonne vis two-nucleon potential [l7| , 
and the Urbana IX three-nucleon potential [l8| : 



ff = E^+E^+ E 



Vi 



(21) 



i<j<k 



A VMC calculation finds an upper bound Ey to an 
cigenencrgy Eq of the Hamiltonian by evaluating the ex- 
pectation value of H in a trial wave function, ^v- 



E x 



>E 



(22) 



Parameters in are varied to minimize Ey, and the 
lowest value is taken as the approximate energy. The 
multidimensional integral is evaluated using standard 
Metropolis Monte Carlo techniques [l9| , hence the VMC 
designation. A good trial function is given by [2(| 



I* 



>=*n 



i<3 



i + u, 



E & 



l*J> , (23) 



where Uij and U^k are non-commuting two- and three- 
body correlation operators induced by the dominant 
parts of Vij and Vijk, respectively; S is a symmetrizer, 
and the Jastrow wave function tyj is 



\*.,)=l[fc(r. lJ )\<S>A(.r;TT z )) 



(24) 



Here the single-particle ^4-body wave function 
<f> a{J t ', TT Z ) is fully antisymmetric and has the to- 
tal spin, parity, and isospin quantum numbers of the 
state of interest, while the product over all pairs of 
the central two-body correlation / c (^j) keeps nuclcons 
apart to avoid the strong short-range repulsion of 
the interaction. The long-range behavior of f c and any 
single-particle radial dependence in $^4 (which, to ensure 
translational invariance, is written using coordinates 
relative to the center of mass of the s-shell core) control 
the finite extent of the nucleus. For p-shell nuclei, 
there are actually three different central pair correlation 
functions f c : f gs , f sp , and f pp , depending on whether 
both particles are in the s-shcll core (ss), both in the 
p-shcll valence regime (pp), or one in each (sp). 



Ui 



(25) 



p=2,6 



where the are the leading spin, isospin, spin-isospin, 
tensor, and tensor-isospin operators in uy. The radial 
shapes of f c (r) and u p (r) arc obtained by numerically 
solving a set of six Schrodingcr-like equations: two single- 
channel for 5=0, T=0 or 1, and two coupled-channel for 
5=1, T=0 or 1, with the latter producing the impor- 
tant tensor correlations [2l[ . These equations contain the 
bare Vij and parametrized Lagrange multipliers to impose 
long-range boundary conditions of exponential decay and 
tensor/central ratios. 

Perturbation theory is used to motivate the three-body 
correlation operator 



Uijk = -eVijk(rij,rjk,rki) 



(26) 



where r = yr, y is a scaling parameter, e is a (small neg- 
ative) strength parameter, and Vijk includes the domi- 
nant short-range repulsion and anticommutator part of 
two-pion exchange in the three-nucleon potential. Con- 
sequently, Uijk has the same spin, isospin, and tensor 
dependence that Vijk contains. 

The variational parameters in / ss , Uij, and Uijk have 
been chosen to minimize the energy of the s-shell nucleus 
4 Hc. For the p-shell nuclei 6 He and 8 Hc, these param- 
eters are kept fixed and the additional parameters that 
enter f sp , f pp , and the single-particle radial behavior of 
3>A have been adjusted to minimize the energy of these 
systems subject to the constraint that the proton and 
neutron rms radii are close to those obtained from more 
sophisticated Green's function Monte Carlo (GFMC) cal- 
culations HHm. 

The wave function samples used here are generated 
by following a random walk guided by the VPy for each 
nucleus. After an initial randomization, a move is at- 
tempted, where each particle is randomly shifted within 
a box of 1.2-1.4 fm in size; the *f?v is evaluated and com- 
pared to the previous configuration, with the move being 
accepted or rejected according to the Metropolis algo- 
rithm. After ten attempted moves, the configuration is 
saved, including the x, y, z coordinates of each particle 
(the center of mass is set to zero) and the probability for 
each particle to be either a neutron or a proton. The 
size of the box gives an acceptance rate of ~50% and we 
generate one million configurations for each nucleus. 

The root-mean-square radii calculated from the VMC 
wave functions are listed in Table HI The second col- 
umn shows the proton root-mean-square radii extracted 
from the charge radii that are obtained from the exper- 
iment based on laser spectroscopy [23j]. The results of 
the VMC wave functions excellently reproduce the radii 
determined experimentally. 

The density distributions of He isotopes are displayed 
in Fig. [TJ Solid and dotted lines denote the neutron and 
proton distributions, respectively. Both proton and neu- 
tron distributions of 4 Hc arc almost the same and arc 
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TABLE I. Root-mean-square radii, given in units of fm, of 
the proton and neutron distributions of 4,6,8 He. The proton 
radius of the second column is obtained by converting the 
measured charge radius. 



He isotope 


Expt. [23] 




Calc. 




proton 


proton 


neutron 


4 He 


1.457(4) 


1.447 


1.447 


6 He 


1.938(23) 


1.928 


2.871 


8 He 


1.885(48) 


1.884 


2.901 



TABLE II. Strength parameters V so of the spin-orbit poten- 
tials for p+ 4 ' 6 ' 8 He scattering at 71 MeV. 



He isotope 


Real potential 


Imaginary potential 


4 He 


0.125 


-0.075 


(i He 


0.025 


0.00 


8 He 


0.05 


0.00 



A. Scattering at 71 MeV 




r(fm) 



FIG. 1. (Color online) Density distributions of 4 ' 6,8 He calcu- 
lated with the VMC wave functions. Solid and dotted lines 
denote the neutron and proton distributions, respectively. 



confined to the small region. Though the proton distri- 
bution is very similar in both 6 He and 8 He, the tail of 
6 He is slightly more extended to larger distances. Be- 
yond 2 fm, the neutron distribution overwhelms that of 
proton. The falloff of 6 He neutron density is rather slow 
compared to that of 8 He owing to its weak binding. 



III. RESULTS 

The construction of the optical potential, © and 
(|13[) . requires the 3A-dimcnsional integration indicated 
in Eq. © or (JTUJ). This integration can be performed 
with ease in the Monte Carlo method using the wave 
function samples generated in Sec. Ill CI The feasibil- 
ity and the accuracy of the Monte Carlo integration in 
the calculation of the phase shift function was already 
demonstrated for various cases in Ref. jHI- The dif- 
ferentiation in Eq. is facilitated by fitting xe(&) in 
terms of several Gaussians with different falloff parame- 
ters: XE(b)=J2k C k^ xl p(~ a kb 2 ). 



First we discuss the elastic scattering observables at 
71 MeV. The strength V^ of the spin-orbit potential is 
the only parameter in the present formalism. The differ- 
ential cross section is largely determined by the central 
potential U c . The spin-orbit potential contributes to the 
differential cross section almost negligibly within a rea- 
sonable range of the strength, and so we try to fit the 
vector analyzing power by varying V^ . In the case of 
4 He, we choose different values of the strength parameter 
depending on the spin-orbit real or imaginary potential. 
The spin-orbit imaginary potential for 6 He and 8 Hc is 
set to be zero for the sake of simplicity. The depth pa- 
rameters of the spin-orbit potential are listed in Table [TT1 
The value of V so for the p+ 4 Hc spin-orbit real potential 
is consistent with that employed in the systematics of 
one-particle motion 25[. Compared to 4 He, the strength 
of the spin-orbit real potential for 6,8 He is considerably 
smaller, in accordance with the analysis of Ref. 

Figure [2] shows the optical potentials calculated for 
p+ 4 ' 6,8 He elastic scattering. The left and right panels 
of Fig. 2(a) are the central real and imaginary poten- 
tials. Solid lines denote potentials of the full calculation, 
while dash-dotted lines are the folding potentials. The 
optical potential for 4 He is much deeper than the others 
at short distances due to its compact structure but be- 
comes much shallower near the surface. Compared to the 
folding potential, the full optical potential has the follow- 
ing property near the surface: the imaginary part turns 
out to be much more absorptive and the real part is less 
attractive. This is a general feature of the DPP due to 
the breakup effect as already observed in Ref. |j[ . Simi- 
larly Fig. 2(b) exhibits the spin-orbit real and imaginary 
potentials. 

Figure [3] displays the differential cross sections (upper 
panels) and vector analyzing powers (lower panels) for 
4 He, 6 He and 8 He calculated using the above potentials. 
Solid lines denote the results based on the potential de- 
rived from the full phase shift function (jTU)) . dotted lines 
are the results obtained from its leading order in the cu- 
mulant expansion (|12[) . and dash-dotted lines are the re- 
sults of the folding potential (fTT| . It is seen that the 
full potential gives smaller and better differential cross 
sections than the folding potential. The vector analyzing 
powers are also better reproduced with the full potential. 
The leading order approximation seems to be surprisingly 
good. The vector analyzing power available for 4 Hc is ac- 
curately known Q, and its behavior that changes a sign 
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(a) Central potentials (b) Spin-orbit potentials 



FIG. 2. (Color online) Optical potentials for p+ 4,6,8 He elastic scattering at 71 MeV. The central and spin-orbit potentials are 
shown in Fig. 2(a) and 2(b), respectively. The left and right panels are the real and imaginary potentials, respectively. The 
spin-orbit imaginary potentials of 6 He and 8 He are set to be zero. Solid and dash-dotted lines are the potentials of the full 
eikonal model and the folding model, respectively. 



around 60 degrees demands a nonvanishing V^ value for 
the imaginary potential. In Ref. @ the effects of a DPP 
due to pickup coupling has been studied on the differen- 
tial cross section and vector analyzing power for proton 
elastic scattering from 6 He. Their results for full calcu- 
lation including both pickup and breakup contributions 
have shown a decrease in the differential cross section 
around 50 degrees and an increase in the vector analyz- 
ing power in angles larger than 50 degrees, which is in 
disagreement with the data. 

Though the present theory reproduces the experimen- 
tal data - especially the vector analyzing powers - fairly 
well, the calculated differential cross sections tend to be 
large compared to experiment. The incident energy of 
71 MeV is likely not to be high enough for the Glauber 
approximation. In such a case that the incident energy 
of the projectile is comparable to its Fermi energy, the 
use of the information on the free space nucleon-nucleon 
interaction may not be valid, but Pauli-blocking effects 
may become important. In fact the chance for the in- 
cident nucleon to collide with the proton target will be 
suppressed because of the Pauli effect. 

To simulate this Pauli-blocking correction phenomeno- 



logically, we replace the a l °^ value of Eq. © with 



-tot „tot f 1 



>pN 



7E F 
5~B 



(27) 



where the Fermi energy Ep is related, in the local den- 
sity approximation, to the nucleon density pn(t) by 
E F = (Sti -2 W( r )) 2/,a '■ The value of 7 is unity ac- 
cording to Ref. [28| . but here it is adjusted to reproduce 
the differential cross section at forward angles. In fact the 
use of Eq. (f27|) with 7 = 1 for 4 He turns out to lead to 
unphysical cross sections because the 4 He density is very 
large and changes drastically in a small region and thus 
the local density approximation may not work so well. 
The adopted value of 7 is 0.5 for all He isotopes. At the 
same time, the range parameter /3 p n is also appropriately 
changed as noted in Ref. [l4[ • The values of V so are kept 
unchanged. The solid line in Fig. [4] denotes the result of 
the calculation with the Pauli-blocking correction. Com- 
pared to the solid line in Fig. [3] that has no correction, 
we see that the Pauli-blocking correction results in the 
decrease of the differential cross sections, leading to fair 
agreement with experiment at forward angles. The vec- 
tor analyzing power for 4 He is very much improved as 
well. Also plotted in Fig. @] are the results of Ref. Q 
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FIG. 3. (Color online) Angular distributions for p+ 4,6,8 He elastic scattering at 71 MeV. The upper panels show the differential 
cross sections and the lower panels the vector analyzing powers. Solid lines are full eikonal calculations (|10[) . dotted lines 
the approximation in the leading-order cumulant expansion H12p , and dash-dotted lines the folding model calculations 
Experimental data are taken from Ref. for 4 He, from Refs. [JH^] for 6 He, and from Refs. @, H3] for 8 He. 
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FIG. 4. (Color online) Angular distributions for p+ 4,6,8 He elastic scattering at 71 MeV compared between calculations and 
experiment. Solid lines are the results of the present model with the Pauli-blocking effect, while dash-dotted lines are the 
results taken from Ref. [f| . See the caption of Fig. [3] for the experimental data. 



that derives the optical potential based on the Watson 
formulation of multiple-scattering theory. Our calcula- 
tion including the breakup effect clearly offers a better 
description of the scattering of the proton with the He 
isotopes. The optical potential with the Pauli-blocking 



correction is displayed in Fig. [5l and it is compared to 
the one without the correction. Both real and imagi- 
nary parts are reduced significantly for r < 2, fm by the 
Pauli-blocking effect. 

As seen in Fig. |U the theory slightly underestimates 
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2 4 2 4 

R(fm) R(fm) 

FIG. 5. (Color online) Optical potentials for p+ 4 He elastic 
scattering at 71 MeV. The upper two panels show the cen- 
tral and spin-orbit real potentials while the lower two panels 
the imaginary potentials. Solid and dotted lines denote the 
potentials of the full calculation with and without the Pauli- 
blocking effect, respectively. 

the p+ 4 Hc differential cross section beyond 70 degrees. 
We considered the effect of knock-on process to see 
whether or not that gives an important contribution in 
improving the cross section. In the knock-on process the 
incident nucleon knocks a nucleon inside the nucleus and 
the knocked nucleon is ejected out of the nucleus. The 
knock-on exchange effect produces a nonlocal potential 
and it was calculated assuming a (0s) 4 wave function 
of 4 He and the Minnesota central potential |2j|. This 
nonlocal potential is transformed to an equivalent local 
potential following the WKB procedure [3fJ. The effect 
of this potential is, however, so small that the above dis- 
crepancy remains to be resolved. 



B. Scattering at 300 and 500 MeV 

Next we study the elastic scattering observables at 
300 and 500 MeV where p+ 4 He data are available for 
both differential cross sections and vector analyzing pow- 
ers [3ll [32} . The eikonal approximation should work bet- 
ter at these high energies. 

Figure [5] displays the optical potentials for proton- 
clastic scattering from 4 He at 300 and 550 MeV. As al- 
ready known it is seen that with increasing projec- 
tile incident energy the depth of the central real potential 
decreases and the central imaginary potential increases. 
The central real potential calculated in the full model at 
550 MeV is found to be extremely shallow, which results 
in a very small real part of the spin-orbit potential if V so 
is of the order of 0.1. In order to account for the vector 
analyzing power, V so has to be chosen to be about 1.0 as 




2 4 2 4 



R(fm) R(fm) 

FIG. 6. (Color online) Optical potentials for p+ 4 He elastic 
scattering at 300 and 550 MeV. The upper two panels show 
the central and spin-orbit real potentials while the lower two 
panels the imaginary potentials. Solid and dash-dotted lines 
denote the full and folding potentials at 300 MeV, whereas 
dotted and dash-dot-dotted lines the full and folding poten- 
tials at 550 MeV. 



will be shown below. 

Figure [7] exhibits observables of proton elastic scatter- 
ing from 4 He at 300 MeV (left side) and 550 MeV (right 
side). From the top, differential cross section (a), vec- 
tor analyzing power (b), and spin-rotation function (c) 
are give n. Dots are experimental data [3l[ for 300 MeV 
and [H for 500 MeV. Note that the theory at 550 MeV 
is compared to the experiment at 500 MeV. The present 
theory reproduces both the differential cross section and 
the vector analyzing power reasonably well. 

Finally we predict the scattering observables at 
300 MeV for p+ 6 Re in Fig. [8] and p+ 8 Re in Fig. M Solid 
and dotted lines in both figures show results for full and 
folding model calculations. In the case where the real 
depth parameter V^o is chosen to be the same as that of 
4 He case at 300 MeV, the angular distributions for 6 He 
and 8 He are similar to that for 4 He at 300 MeV, as ex- 
pected. Calculations with two times larger parameter 
(V so = 0.2) increase the vector analyzing powers at small 
angles and the differential cross sections at large angles. 
The differential cross section, the vector analyzing power, 
and the spin-rotation function arc all different between 
the full and folding model calculations at angles larger 
than 50 degrees, as we observe in Figs. [3] and [7] 



IV. CONCLUSION 

We have analyzed the elastic scattering observables 
for protons scattered from the He isotopes at 71 and 
300 MeV. The optical potentials for p+ 4,6,8 He systems 
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FIG. 7. (Color online) Angular distributions for p+ 4 He elastic 
scattering at 300 and 550 MeV. Displayed are the differential 
cross section (a), the vector analyzing power (b), and the spin- 
rotation function (c). The values of V so are 0.1 and —0.05 for 
the real and imaginary potentials at 300 MeV, and 0.9 (0.09 
for folding calculation) and —0.025 at 550 MeV. Solid, dotted, 
and dash-dotted lines indicate the same types of calculations 
as those of Fig. [3] Experimental data are taken from Ref. [3l[ 
for 300 MeV and from Ref. ^ for 500 MeV. 



FIG. 8. (Color online) Angular distributions for p+ 6 He elastic 
scattering at 300 MeV. The left panels show the results calcu- 
lated with 1/50=0.1 for the real potential, the same strength 
as that of 4 He at 300 MeV, while the right panels the results 
with V^ o =0.2. The spin-orbit imaginary potential is set to be 
0.05, which is the same as that of 4 He. Solid, dotted, and 
dash-dotted lines indicate the same types of calculations as 
those of Fig. 



are calculated in the Glauber model. The central poten- 
tial is evaluated to all orders of the complete Glauber am- 
plitude using the nucleon-nucleon scattering amplitude 
and the ground-state wave function of the He isotope 
that is taken from the variational Monte Carlo method. 
Both the real and imaginary parts of the central poten- 
tial are determined without any adjustable parameters. 
It should be noted that the central potential obtained in 
this way takes into account the breakup effect of the He 
isotope to its excited states including continuum states, 
which makes it possible to learn the difference from the 
single folding model potential. The spin-orbit potential is 
assumed to take the standard form that uses the deriva- 
tive of the central potential. Its strength is the only pa- 
rameter in the present approach. 

Though the incident energy of 71 MeV may be a lit- 
tle too low for applying the Glauber theory, the present 
theory leads us to reasonable agreement with experimen- 



tal data especially on the vector analyzing powers for 
p+ 4 ' 6,8 Hc scatterings simultaneously. It should be noted 
here that usual i-folding calculations fail to reproduce the 
vector analyzing powers. We observe that the differential 
cross sections are all slightly larger than the experimen- 
tal data even at forward angles. We have studied the 
Pauli blocking effect that partly suppresses the interac- 
tion between the proton and the He isotope and found 
that including its effect improves the angular distribution 
of the elastic scattering. 

For higher incident energy scattering from 4 He, the 
differential cross section is reproduced well at forward 
angles and the vector analyzing power is also reproduced 
well in accordance with the experimental data. The vec- 
tor analyzing power in the intermediate energies has been 
a long standing problem that defies reproduction of the 
experimental data. However, the spin-orbit potential 
given by the derivative of the central potential calculated 
from the Glauber theory appears to be able to reproduce 
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the proton-elastic scattering in both low and intermedi- 
ate energy region. We have predicted the angular distri- 
butions for p+ 6,8 He elastic scattering at 300 MeV. 



FIG. 9. The same as Fig. [8] but for p+ 8 He elastic scattering 
at 300 MeV. 
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